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The vacuum of Quantum Chromodynamics is a superposition of degenerate states with different 
topological numbers that are connected by tunneling (the 6*- vacuum). The tunnehng events are 
due to topologically non-trivial configurations of gauge fields (e.g. the instantons) that induce local 
"P-odd domains in Minkowski space-time. We study the quark fragmentation in this topologically 
non-trivial QCD background. We find that even though QCD globally conserves V and CV symme- 
tries, two new kinds of P-odd fragmentation functions emerge. They generate interesting dihadron 
correlations: one is the azimuthal angle correlation ~ cos((/>i +412) usually referred to as the Collins 
effect, and the other is the P-odd correlation ~ sin(0i -I- (^2) that vanishes in the cross section 
summed over many events, but survives on the event-by-event basis. Using the chiral quark model 
we estimate the magnitude of these new fragmentation functions. We study their experimental 
manifestations in dihadron production in e^e~ collisions, and comment on the applicability of our 
approach in deep-inelastic scattering, proton-proton and heavy ion collisions. 



1. Introduction. Quantum Chromodynamics (QCD) is 
at present firmly established as the theory of the strong 
interactions. Equations of motion in QCD possess topo- 
logically non-trivial solutions [l| signaling the presence of 
degenerate ground states differing by the value of topo- 
logical charge 0. The physical vacuum state of the the- 
ory is a superposition of these degenerate states, so-called 
^-vacuum Q. To reflect this vacuum structure one may 
equivalently introduce a 0-term in the QCD Lagrangian. 
Unless 9 is identically equal to zero, this term explicitly 
breaks V and CV symmetries of QCD. However stringent 
limits on the value of < 3 x 10"^*^ deduced from the 
experimental bounds on the electric dipole moment of 
the neutron Q indicate the absence of global V and CV 
violation in QCD. 

Nevertheless it has been proposed that the local V- and 
CP-odd effects due to the topological fluctuations char- 
acterized by an effective 9 = 9{x, t) varying in space and 
time could be directly observed through multi-particle 
correlations In heavy ion collisions, the existence of 
magnetic field (and/or the angular momentum) in the 
presence of topological fluctuations can induce the sep- 
aration of electric charge with respect to the reaction 
plane, so-called Chiral Magnetic Effect ^-^\. There is a 
recent experimental evidence for this effect from STAR 
Collaboration at RHIC [lo| . The interpretation of STAR 
result in terms of the local parity violation is under in- 
tense scrutiny at present, see e.g. 11- isj. 

In this paper, we study the role of QCD topology in 
hard processes using the formalism based on factorization 
theorems IJ]. From the QCD factorization point of view, 
the cross section in high energy collision can be factor- 
ized into a convolution of perturbatively calculable par- 
tonic cross section and the non-perturbative but univer- 
sal parton distribution and fragmentation functions. In 
the conventional formalism, these distribution and frag- 
mentation functions are required to be P-even because 



of the parity-conserving nature of the strong interaction. 
However, in the presence of local (in space and time) P- 
odd domains P-odd fragmentation functions can emerge 
15j ; note that only the cross section of the entire process 
has to be P-even, not the fragmentation function. 

In this letter, we derive the most general form of the 
quark fragmentation function for a quark fragmenting 
into a pseudoscalar meson which is consistent with the 
Lorentz invariance. Abandoning the parity constraint, 
we obtain two P-odd fragmentation functions besides the 
well-known P-even spin-averaged fragmentation function 
[lit ^^'^ Collins function (T^]. We obtain the exact op- 
erator definitions and estimate the size of these new P- 
odd fragmentation functions using the chiral quark model 
As a first step, we present their observable effect in 
the back-to-back dihadron production in e+e^ collisions. 
We encourage the experimentalists to carry out the re- 
lated analyses at RHIC and elsewhere. 

2. Quark fragmentation functions in locally V-odd 
background. The quark fragmentation functions are de- 
fined through the following matrix [l9| : 



1 r dy d^yjL ,k.y 
z J (2^)3 

{PX\4,{0)4my^=0, 



Cyi^{v)\PX) 



(1) 



where p is the momentum of the final state hadron with 
a transverse momentum p±^ relative to the fragment- 
ing quark k. We choose the hadron moving along -\-z 
direction, and define the light-cone momentum = 
(/ ± p^)/V2. For convenience, we define two light- 
like vectors: = 5^+ and = . The momen- 
tum fraction z — p'^ /k^, and k± = —p±/z. Cy = 
V exp [ig dXn ■ A{y + An)) is the gauge link needed 
to make A{z,p±) gauge invariant. 

Since QCD is a theory conserving C, P, and T globally, 
one usually expands the above matrix using the following 
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constraints [20|: 

Hcrmiticity : 
Parity : 
Time — reversal : 



At(p,fc)=70A(p,fc)7" (2) 
A(p,fc)=70A(p,fc)70 (3) 
A*(p,fc) = FTA(p,fc)Ffi(4) 



where Vr 



ij^j^ and 



(p", — p). Using the 



basis of gamma matrices T — {1,7^, 7^7^, ct^", 17^}, and 
the available momenta p and k, one can expand A(p, k) 
in the most general form: 



Aip.k) 



M A1I + + A3ft + A^a 



M 



Ar^i)-I^ + Af,^-i^ + MAji-i^ 



(5) 



where Af is the hadron mass used to make all A^'s have 
the same dimension. Since the time-reversal changes out- 
state to in-state, it does not really give any constraint 
on the coefhcients Ai 20]. One the other hand, if one 
applies the Hermiticity constraint, all of the A^'s have 
to be real. If one further applies Parity constraint, one 
finds A5 — Aq = Ar = As = 0. However, as we stated 
in the Introduction, we are interested in the situation in 
which a local 7^-odd domain develops in space-time, and 
the quark fragmentation happens inside such a T'-odd 
domain (or in other words, the quark scatters off the non- 
trivial gauge field configuration prior to transforming into 
a pseudoscalar meson). In this case, the T'-odd modes 
in the quark fragmentation could be populated and one 
has to release the parity constraint in Eq. ([3]). Without 
parity constraint, we thus need to keep all 8 terms Ai 
through As in Eq. ([5]). Applying the twist-expansion 
by parametrizing the momenta as w p^n^ and k^^ w 



p'^)/z and keeping the leading terms we obtain 



Mz,P±) 



1 r 
2 



D{z,pl)f. + HHz,pl)a^'^P-^ 



+ - D{z,pIW 
+Ht{z,p\)a^^''i^'- 



M 



(6) 



where D{z,p\) and H^{z,p']_) are the usual T'-even 
fragmentation functions: D{z,p\) is the transverse mo- 
mentum dependent spin-averaged fragmentation function 
0, and Ht{z,p\) is the Collins function describing 
a transversely polarized quark fragmenting into an un- 
polarized hadron (17f. Now besides the two conven- 
tional 7^-even fragmentation functions, we also obtain 
two new T'-odd fragmentation functions: D{z,p\) and 
Hi{z^p\). As we will show below, H^{z,p\) generates 
a new kind of azimuthal correlation. Its role is simi- 
lar to Hi{z,p\): H^{z,p'^) represents an asymmetric 



distribution oc {p x p±) ■ Sq, while H^{z,p'\^ represents 
an asymmetric distribution oc p± ■ Sq for a transversely 
polarized quark with spin vector Sq to fragment into a 
pseudoscalar meson. The newly derived 7^-odd fragmen- 
tation functions will lead to interesting P-odd effects in 
experiment as we will show in the next section. 

In order to study the experimental effects generated 
by these P-odd fragmentation functions, we need to es- 
timate their magnitude. For this purpose we use the 
effective chiral quark model developed by Manohar and 
Georgi which is an effective theory of QCD at low 
energy scale. This model has also been adopted for an 
estimate of the Collins functions in [2l|, [l^] . The effec- 
tive Lagrangian describing the interaction between the 
quarks and the pion in the leading order is given by 

Lggn = - If V5?7^7 V • d^TTi^q (7) 
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where /^r ~ 93 MeV is the pseudoscalar decay constant. 



FIG. 1: Lowest-order Feynman diagram for a quark with mo- 
mentum k fragmenting into a n meson with momentum p. 

At tree level, the fragmentation of a quark is modeled 
through the process q* — > irq, see Fig. [T] One can obtain 
the unpolarized quark fragmentation function D{z,r)]) 
from the definition in Eq. ([6|) which has been done in [2l| . 
The Collins function can be calculated similarly, though 
one needs to go beyond the tree diagram and consider 
the TT-loop to obtain the final result, see Ref. 
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Since the chiral quark model Lagarangian in Eq. ([7]) 
conserves parity, it does not generate 7^-odd fragmenta- 
tion functions D{z,p\) and Hi{z,p\). As we stated 
in the Introduction, QCD contains topological gauge 
field configurations, and their effect can be mimicked by 
an effective space-time dependent 6 field [1, @]. One 
can thus add to the Lagrangian of QCD the term 
{g'^ /'i2TT'^)e{x,t)F^'' F^^; performing an axial C/(l) rota- 
tion this term can be transformed into 2W]^i^^'^q'y^^^'^q 
Q. Let us define an effective 9^ — d^d/2Nf, whose zero 
(time) component is the chiral chemical potential /is in- 
troduced in [9] . The existence of this new term will yield 
a modified quark propagator iS{p, 9) = i/(J) — m + ^7^) 
given by 

iS{p,e) = i[VtzS{p + 0)+VcS{p-0)] 

X [i + m-f'^ {s{p + e)-s{p-e))] 



iip + oy 



-1 
(8) 
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where Vc.Ti are the left (right) projection operator 
Vc,n = (1 ± 7^)/2, and iS{p) = i{i> + m)/{p'' - m^) 
is the conventional quark propagator. 

With the modified quark propagator in the 7^-odd 
"bubble" , both the T'-odd fragmentation functions can be 
generated at the tree level as shown in Fig.[TJ We further 
find that the new fragmentation functions are suppressed 
by a factor ^""^/p"*" when 9 is along time or z-direction. 
Since p"*" is a large component in our twist-expansion, we 
thus neglect the effect from the 0, 3 components of 6 to 
be self-consistent. On the other hand, if the 9 is along 
the transverse direction that is perpendicular to p±, we 
find that the T'-odd fragmentation functions vanish. We 
thus only consider the situation when 9 is along the p± 
direction: 0^ — 9±p^, in which case we find 



D{z,pI) 



9\ 



4:9±P_L 



1 - - 



(9) 



hHz,pI) 



GAf^TT^z p'^ + z'^m-^ + (1 — z)m'^ 
4(1 — z)^z^m^ 
(p^ -I- z^m'^ + (1 — 

g\ mqm^ 9j_ 1 

4/2 87r3 P±(pl + z^ml + {l-z)miy 

{pl + z^mlf{z^2) + {l-zfml 



[(3z - 2)r 



4(p5 



2 2 

z m 



(10) 



where we have neglected terms ~ 0{9^). From above 
equations, we see that both of the T'-odd fragmentation 
functions are proportional to 9±/pj_. Since p± is a small 
component, the effect is not suppressed. We will now 
estimate the size of the observable effect generated by the 
T'-odd fragmentation functions within the same model. 

3. Observable effect of parity-odd fragmentation func- 
tions. Let us now discuss the experimental consequences 
of the V-odd fragmentation functions. As a first step, we 
study a relatively simple process, the back-to-back di- 
hadron production in e+e" collisions e+e~ — >■ + X. 
The method we presented here can be generalized to 
study T'-odd effects in heavy ion collisions. 

At leading order in QCD coupling, the two hadrons hi 
and /i2 in e~^e~ collisions are the fragmentation products 
of a quark and an antiquark originating from e^e~ — qq 
annihilation. Following Ref. [23|, we choose a reference 
frame such that the e"^e~ — >■ qq annihilation occurs in the 
x-z plane, with the back-to-back quark and antiquark 
moving along the z-axis. The final hadrons hi and /12 
carry light-cone momentum fractions zi and Z2 and have 
intrinsic transverse momenta pi± and p2± with respect 
to the directions of the fragmenting quarks. Using the 
fragmentation parameterization in Eq. ([5]) , one can derive 



the differential cross section as 
da 



dVS 



'^o5Ie^{(l + cos'0) 
X \Dq{zi)Dq[z2) - Dg{zi)bq{z2) 



- sin 9 cos(©i -I- ( 



X H^{zi)H^{z2) + H^{zi)H^{z2) 
+ s\v? 9 sin((/)i -I- 02 ) 

X 'h^{zi)H^{z2)-H^{zi)H^{z2) 



:ii) 



where the phase space dVS = dzidz2d cos 9 d{4>i -\- 4'2), 
Co — NcCt1jn/4:Q'^, and 9 is the angle between the ini- 
tial beam direction and the z-axis, not to be confused 
with the 9{x) field. In Eq. (|lip . we have integrated over 
the moduli of the intrinsic momenta, pi± and P2±, and 
over the azimuthal angle (j)i . The p^-integrated functions 
Dq{z) and H^{z) are defined as 



d'p^Dqiz,pl), (12) 
dW^MHt\z,pl). (13) 



The definition of Dq{z) (or H^{z)) is similar to Dq{z) 
(or H;^iz)). 

The cos((/)i -I- 02) correlation is usually referred to as 
the Collins effect, analyzed recently by BELLE CoUab- 
However, we find that the product of two 
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oration 

7^-odd fragmentation functions Hq{z) leads to the same 
azimuthal correlation. This complicates the extraction 
of the Collins function, but may in effect provide an al- 
ternative view of the origin of the Collins effect and puts 
an experimental constraint of the V-odd fragmentation 
function H^{z). It is interesting that a new azimuthal 
correlation also emerges: the sin((/)i -1- 02 ) term, which 
is explicitly V-odd. Notice that for the sin(0i -1-02) 
contribution, the first term corresponds to the situation 
when the antiquark fragments inside the P-odd bubble, 
see Fig. HJa), whereas the second term corresponds to 
the situation when the quark fragments inside the T'-odd 
bubble, see Fig. [^b). They have the opposite sign, and 
thus when averaged over the large number of events, the 
effect will vanish. Thus a 7^-odd effect happens only on 
the event-by-event basis Q. 



(a) 



(b) 



FIG. 2; Illustration of the antiquark (a) and the quark (b) 
fragmenting in the P-odd bubble. 
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To estimate the effect, let us assume that the anti- 
quark fragments inside the 'P-odd bubble; the relative 
magnitude of the correlation will depend on the fol- 
lowing factor 1(9, zi, Z2), besides the kinematic factor 
sin^ 9/{l + cos^ 6*), 

I{0, zu Z2) = ^ ^ ^ V^' ^ . (14) 

Dq{z^)Dq{z2)-Dq{zi)Dr,{z2) 

Certainly I{9, zi, 2:2) depends on the size of 9^. To esti- 
mate 9±, we resort to the instanton vacuum model (for 
a review, see [25]). According to [2^, the two most im- 
portant parameters are the mean size of the instanton 
p '--^ 1/3 fm and the typical separation R between instan- 
tons which satisfies p/R ~ 1/3. We thus estimate 9±_ as 
follows: 



1 



2N 



f 



■d±9{x,t) 



1 1 
2Nf'p' i?2 



10 MeV, (15) 



where we have used Nf = 3, and the factor p^/R^ repre- 
sents the probability for a quark scatters off an instanton; 
note that in Minkowski space-time the instanton event is 
elongated along the light cone [l3|. With = 10 MeV, 
and other standard parameters of the chiral quark model 
[l8| , and using the calculation of the fragmentation func- 
tions taken from Ref. [2l|, we find I{9,zi,Z2) ~ 1.5% for 
a typical zi — Z2 = 0.5 at BELLE experiment, with the 
final two hadrons as 7r+ and 7r~. We urge the experi- 
mentalists at BELLE, RHIC and elsewhere to carry out 
an analysis to constrain the P-odd fragmentation func- 
tions. Because of the universality of the fragmentation 
functions, we expect that the formalism developed here 
could be generalized to other processes. 

4- Conclusion. In this letter we have studied the quark 
fragmentation in the topologically non-trivial QCD back- 
ground. We have found two new fragmentation func- 
tions besides the well-known spin-averaged fragmenta- 
tion function and the Collins function. Both of the new 
fragmentation functions are P-odd. We have related the 
magnitude of these functions to the typical size of the 
topological fluctuations (described by the effective 9{x) 
field) . We have studied the observable effects of the V- 
odd fragmentation functions in back-to-back dihadron 
production in e+e^ collisions, and have found that a new 
azimuthal correlation oc sin((/)i -I- </)2) appears. Since the 
new azimuthal correlation is explicitly P-odd, it can be 
observed only on an event-by-event basis. Our results 
also offer a new interpretation of the Collins correlation. 
We encourage the experimentalists to carry out an anal- 
ysis to constrain the P-odd fragmentation functions, and 
anticipate new applications, e.g. in heavy ion collisions. 
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